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The magnetic properties of single-domain nanoparticles with different geometric shapes, crys- 
talline anisotropics and lattice structures are investigated. A recently proposed scaling approach is 
shown to be universal and in agreement with dimensional analysis coupled with an assumption of 
incomplete self-similarity. It is used to obtain phase diagrams of magnetic nanoparticles featuring 
three competing configurations: in-plane and out-of-plane ferromagnetism and vortex formation. 
The influence of the vortex core on the scaling behavior and phase diagram is analyzed. Three- 
dimensional phase diagrams are obtained for cylindrical nanorings, depending on their height, outer 
and inner radius. The triple points in these phase diagrams are shown to be in linear relationship 
with the inner radius of the ring. EUiptically shaped magnetic nanoparticles are also studied. A 
new parametrization for double vortex configurations is proposed, and regions in the phase diagram 
are identified where the double vortex is a stable ground state. 



I. INTRODUCTION 

Magnetic thin fihus and nanoparticles have been in- 
tensively studied during the last two decadesii^, not only 
because of their great potential for technological appli- 
cations, but also because of fundamental scientific inter- 
est. Many new phenomena come about by imposing ge- 
ometric restrictions in onei or more dimensions^ii^i^i^. It 
has been demonstrated that nanoparticles show predom- 
inantly single domain structure when their size is smaller 
than a characteristic length scaled. These single do- 
main particles are promising candidates for high density 
data storage^., integrated magnetic-electronic devices^, 
and applications in biotechnologyi^. 

A great deal of attention has focused on arrays of 
magnetic nanoparticles. In magnetic nanoparticle ar- 
rays, there are two distinct issues of interest: the spin 
configuration of the individual particles, and the inter- 
actions between them. Here we focus on the first issue. 
The magnetic properties obtained under this considera- 
tion are valid when the distance between the individual 
particles is larger than twice the characteristic size of the 
individual particles since it has been shown that the inter- 
actions can be safely neglected under this conditioui^iii 

Within nanoparticles, different magnetic configura- 
tions have been observed, including vortex, leaf, and 
flower statesji^ Single-domain configurations have at- 
tracted continuous attention for their obvious application 
potential. In particular, the magnetic vortex, also known 
as non-localized soliton has been explored recently for its 
application potential and interesting dynamics.— li^*^ In 
this work, wc study magnetic phase diagrams of such 
nanostructures as a function of their shape, crystalline 
anisotropy and lattice structure. 

On the numerical side, a scaling approach has been 
shown to be effective in determining phase diagrams for 
cylindei— and cone>il shaped nanoparticles. Here, we 
provide a systematic numerical study for different ge- 
ometric shapes. When the characteristic length scale 
is sufficiently small, the shape of the particle is one 



of the dominant factors determining its magnetic prop- 
erties. Numerous experimental investigations have ad- 
dressed this and related issues of domain structurej^i^ 

On the conceptual side, the scaling approach s ugg ests 
a self-similarity of magnetic nanoparticles Ref. [20 tac- 
itly assumes that magnetic nanoparticles exhibit com- 
plete self-similarity with respect to the small parame- 
ter a/ Lf^x, where, a is the lattice spacing and L(,x is the 
magnetic exchange length. In fact, the particles exhibit 
only incomplete self-similarity with respect to the lattice 
spacing in certain circumstance, as we demonstrate in 
this work. This incomplete similarity agrees with dimen- 
sional analysis, as it must, and with available numerical 
data^^ 

The topology of the nanoparticle plays an important 
rolei^i^i^i^ In a simply connected topology, vortex 
states, for example, typically must have a core region in 
which the spins point out of the vortex plane. In a nanor- 
ing, furthermore, the inner radius Ri provides an addi- 
tional length-scale with which to probe the self-similarity 
of magnetic nanoparticles. We perform a scaling analy- 
sis, and show that nanoparticles in this topology exhibit 
complete similarity with respect to the lattice constant. 
This is a consequence of the additional length Ri that 
plays the role of the lattice constant in regulating the 
vortex core energy. 

Quite a different scenario of the magnetization rever- 
sal was revealed in elliptical particles . Numerical 
simulations showed different spin configurations includ- 
ing multi vortex states^. The double vortex configura- 
tion confined in cUiptically shaped ferromagnetic particle 
is especially interesting, for it provides a model system for 
the study of static and dynamic interaction between soli- 
tons (localized solution of nonlinear equations)^. Many 
efforts^Si^ have been taken to obtain the phase diagram 
of such systems, yet it is still an open problem. 

The first topic of this paper is to reveal the essence of 
the scaling approach and to verify its validity in terms 
of different shapes, anisotropy and crystalline structure. 
Also, the effects of these parameters on the phase dia- 
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grams are analyzed. The influence of the vortex core on 
the scahng behavior and phase diagram is investigated. 
Furthermore, the scahng approach is apphed to nanor- 
ings and eUiptically shaped nanoparticles. The resulting 
phase diagrams are given, and new and interesting phe- 
nomena are discussed. 



II. MODEL AND NUMERICAL PROCEDURE 

In the absence of an external magnetic field, the Hamil- 
tonian (7i)(or energy) of a magnetic nanoparticle consists 
of three terms: exchange interaction, dipolar interaction, 
and crystalline anisotropy. If each magnetic moment oc- 
cupies a site of the underlying lattice, Ti. is given by 

<i,j> 

+dY, ^'' ~ ^^^^3 ' ^''^ + Uk, (1) 

where J > is the ferromagnetic exchange constant (or 
exchange integral, measured in units of energy)^, which 
is assumed to be non-zero only for nearest neighbors (nn) , 
D is the dipolar coupling parameter and fij the displace- 
ment vector between sites i and j. The anisotropy term 
Uk can take various formal among which the most com- 
mon arc uniaxial anisotropy Uk — K^^sin^Oi , where 
9i is the angle Si makes with the easy axis, and cu- 
bic anisotropy Uk = K^,^[af0f + Pfjf + aj^f], where 

<^i,Pi,Ji are the direction cosines of Si. Note that K 
is the single site anisotropy energy (not an energy den- 
sity). For most materials, the dimensionless ratio D/ Ja? 
falls in the range of 10"'^ and 10^*, where, the lattice 
constant a is approximately SA. The dimensionless ratio 
Ka^ /D lies between and 10. We choose Ka^ /D — 1, 
Ja^ /D = 5000 and a = 3 A in the following calculations, 
unless they are specified otherwise. 

The objects studied in this paper are magnetic 
nanoparticles with various shapes and anisotropics. In 
such systems, three dominant competing configurations 
have been identifiedi^: (I) out-of-plane ferromagnetism 
with the magnetization aligned parallel to the nanodot 
base; (II) in-plane ferromagnetism with the magnetiza- 
tion perpendicular to the base; (III) a vortex state with 
the magnetic moments circling in the base plane. Dou- 
ble vortex states in eUiptically shaped particles will be 
discussed in detail later. A typical phase diagram for a 
cylinder is shown in Fig. 1, which exhibits these three 
phases function of the cylinder radius R and its 

height H. Note that there can be other metastable con- 
figurations, such as the buckle state^i^i^ which are not 
considered here. These states result from the competi- 
tion between the exchange and dipolar interaction. The 
exchange interaction tends to align spins in the same di- 
rection, whereas the dipolar interaction encourages spins 



to minimize their magnetostatic energy resulting in the 
shape anisotropy. Thus spins align in-plane in a flat disk 
while they point out-of-plane in an elongated cylinder. 
The vortex state is also a result of dipolar interactions 
since it nearly eliminates the demagnetization field. 

In order to obtain phase diagrams such as the one 
shown in Fig. 1, one could resort to analytical calculation 
based on a continuum model. However this approach is 
limited to highly symmetric shapes and magnetization 
configurations. An alternative is to use numerical sim- 
ulations. These can be powerful and universal but are 
often limited by computational resource. As outlined in 
Ref. [l6l . the major technical problem is that the number 
of magnetic moments in systems of physical interest is of 
the order of 10^, which cannot presently be handled, even 
by high-end supercomputer facilities. To overcome this 
restriction, a scaling approach was recently proposed and 
demonstrated for cylinder— and coneii shaped nanopar- 
ticles. They showed that the phase diagram for an ar- 
tificial small J' = xJ {x < 1) could be scaled to the 
phase diagram for the original J according to L' = x^L 
[rj ~ 0.55 and L can be i?, H). The phase boundary for 
small J' appears at small sizes which involve less number 
of spins so that lots of computing time is saved. 
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FIG. 1: (color online) Scaled phase diagram of cylindrical 
magnetic nanoparticle {Ka'^ /D = 1 and Ja? /D — 5000) as 
a function of its radius and height. The underlying lattice 
is simple cubic with a — 0.3A. The three competing phases 
are (I) out-of-plane ferromagnetism, (II) in-plane ferromag- 
netism, and (III) the vortex state. The transition lines are 
obtained according to the scaling approach discussed in the 
text. The inset shows the dependence of the radius at the 
triple point on the scaling factor x. 

This proposal is equivalent to dimensional analysis 
coupled with a statement of incomplete similarity. We 
seek, for example, to find the height H separating the 
vortex phase from the ferromagnetic phase(s). In addi- 
tion to the two governing parameters J, D appearing ex- 
plicitly in Eq.[Tl we also have the radius, R of the cylinder 
and the lattice constant a. Thus we seek a physical law 
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for the critical height of the following form: 

H^f{J,D,R,a). (2) 

From dimensional analysis, only two of the four govern- 
ing parameters have independent dimensions. Following 
convention, we choose the independent parameters to be 
J and D, define the exchange length L^^. = a^/joF/D, 
and express the scaling law in dimensionlcss form: 

n = $(ni,n2), (a) 

where 11 = H/Lex, Hi = R/ Lf^x,^2 — a-lLexi and the 
scaling function $ does not depend on the governing pa- 
rameters of independent dimension J, D. 

Now the typical values a ~ O.Snm, Lex ^ 20nm give 
112 <C 1. We are tempted to suggest complete similarity 
with respect to the small, dimensionless governing pa- 
rameter 1121^ Hence we consider the limit 112 = 0: 

n = $(ni,o) = $i(ni), (4) 

where $i is independent of J and D and a. Recasting 
in original variables, we have that H = L^x^iiR/ L^x)- 
Notice now the invariance of this relation under the fol- 
lowing rcscaling of the governing parameters and critical 
height: 

J' = xJ, 
D' = D 
R' = x'/^R 

a = a 

H' = x'/^H, (5) 

where, x is any positive number. One way to see this 
is to notice that all lengths entering Eq. [4] get rescaled 
by the same amount x^/^, thus the dimensionless ratios 
n, Hi, arc invariant. The invariance of Eq. [blunder this 
transformation is a consequence of dimensional analysis 
combined with complete similarity with respect to the 
dimensionless governing parameter 112 = ol/ Lex- 

Interestingly, the numerical calculations under the as- 
sumption of a core-free vortex phase, where only the mag- 
netic moment located exactly at the center of the vortex 
has a component pointing out of the vortex plane, do not 
obey this scaling.— Instead, they exhibit only incom- 
plete similarity with respect to 112 = a/^e^'^ Namely, 
for small values of 112 ^ 1 , we have 

n = nr''$2(-i^^), (6) 

where, the constant r/ k, 0.55 does not follow from dimen- 
sional analysis. Here, $2 is independent of J, Z), and a. 
The special case of complete similarity is obtained when 
?/ = 1/2. This implies that the physical law Eq. [2l is not 
invariant under the transformation in Eq. [51 Rather, we 
find invariance under the modified transformation: 



J' = xJ, 

D' = D 

R' ^ x''R 

a = a 

H' = ;r"i/. (7) 

This is precisely the transformation described in Ref. 
[igI . and is a consequence of dimensional analysis com- 
bined with incomplete similarity with respect to the small 
dimensionless parameter n2 = a/ Lex- This incomplete 
similarity results from the fact that there is a singularity 
in the magnetization function. 

Whenever we are presented with a scaling phenomenon 
such as Eq. [SI we have the opportunity to save consider- 
able computational and experimental effort. The scaling 
law expresses a physical similarity between systems with 
different values of the governing parameters, so that we 
can use one to study the other. In particular, the authors 
of Ref. [1^ suggest that we study small systems with small 
exchange constant J, which arc less computationally in- 
tensive to simulate. Then use Eq. [Tjto scale up the results 
to the large systems with large exchange constant that 
are of immediate physical and technological interest. The 
proposal is justified by the incomplete similarity of the 
physical law Eq. [^jwith respect to the small, dimension- 
less parameter H2 = a/ Lex- Other physical quantities of 
nanomagnets may satisfy incomplete similarity, including 
dynamic and thermal properties j^i^^i^. 

III. RESULTS AND DISCUSSION 
A. Shape, anisotropy, and lattice structure 

Incomplete similarity with respect to the lattice con- 
stant occurs in magnetic nanoparticles regardless of 
cross-sectional geometry, crystalline anisotropy, or lat- 
tice structure. To illustrate the use of the scaling pro- 
cedure, let us first consider the example of a cylindrical 
nanoparticle. Using the 2000-nodc 15.78 tcraflop high- 
performance supercomputer at University of Southern 
California (USC), the energies of the competing phases 
were evaluated throughout the parameter plane spanned 
by the cylinder radius R and height H for systems with 
up to 400,000 sites. The scaling procedure was then used 
to collapse the resulting phase diagrams with different 
scaling factors, four of them (x=0.02; 0.04, 0.06, and 
0.08) given in Fig. 1 as examples. Note that there is 
a triple point Ht), which is used to extract the scal- 
ing exponent, shown in the inset of Fig. 1. For the sake 
of simplicity, a simple cubic underlying lattice structure 
with cubic crystalline anisotropy and the "core-free" vor- 
tex state is adopted. Discussion about other structures 
and the effect of the core will come later. 

The scaling exponent 77 = 0.556 is consistent with the 
previous resulti^, suggesting incomplete similarity with 
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respect to the lattice constant in this case. It is observed 
that the slope of the line separating the two ferromag- 
netic phases is k=1.811, which is in exact agreement 
with the analytical solution given previouslj*^ and ar- 
gued later—. 

Since an enormously wide range of magnetic properties 
can be obtained by using different geometric shapes^, it 
is of great interest to see whether nanoparticles with dif- 
ferent cross-sectional geometry exhibit incomplete simi- 
larity as well. To answer this question, here we consider 
prism shaped nanoparticles with triangular, square, pen- 
tagonal, and hexagonal cross sections. From the results 
shown in Fig. 2(a) we find that within an error bar of 
2%, these different geometries have the same scaling ex- 
ponent showing incomplete similarity. In spite of the 
apparently universal scaling behavior, it is also evident 
that different geometries do favor different spin config- 
urations. More precisely, the more symmetric the cross 
section is, the more the vortex phase is favored. Obvi- 
ously, cylindrical nanodot favors the vortex configuration 
the most. Another property of interest is the slope k of 
the line separating the two ferromagnetic (FM) phases. 
Fig. 2(b) shows this slope as a function of the cross sec- 
tion area. To compare the various polygon shapes, they 
have been normalized such that the distance from the 
corner of each polygon to its center is unity. The slope 
is found to increase with the basal area. This trend is 
easy to understand, since the two FM configurations are 
determined by dipolar interactions, i.e. via the demag- 
netizing field which in turn is related to the surface area. 
Quantitatively the slope is expected to be approximately 
proportional to the square root of the area, which is found 
to be in agreement with the numerical results shown in 
Fig. 2(b). 

In the following analysis of the universality of scal- 
ing for various crystalline anisotropy and underlying lat- 
tice structures, we will focus on cylindrical shapes for 
the simple reason that these are most commonly found 
in the existing experimental literature. Fig. 2(c) gives 
phase diagrams for different anisotropics. In accordance 
with intuition, cubic anisotropy favors the two ferromag- 
netic phases equally, i.e. the slope separating these two 
phases does not depend on Ka'^/D, and at the same time 
suppresses vortex formation. Hence, one should consider 
materials with small cubic anisotropy if one wishes to sta- 
bilize the vortex state. Besides cubic anisotropy, another 
prevalent type is the uniaxial anisotropy. This anisotropy 
typically exists in hexagonal close-packed (hep) lattices, 
but it can also occur in cubic lattices due to coupling to 
the substrate or other parts of the environment. In our 
calculation, the easy axis is set to be along the axis of the 
cylinder. The resulting phase diagram is shown in Fig. 
2(c). We observe that uniaxial anisotropy does not affect 
the scaling behavior and exponent. However, a feature 
worth mentioning is that uniaxial anisotropy does change 
the slope of the line separating the two ferromagnetic 
phases, favoring out-of-planc alignment (phase I). The 
larger the value of Ka^/D, the smaller the slope (see Fig. 
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FIG. 2: (color online) (a) Scaled phase diagrams for prism 
shaped nanoparticles. The radii R are defined as the distance 
from the base center to the corner of the polygons. The ex- 
tracted scaling exponents for the triangle (T), the square (S), 
the pentagon (P) and the hexagon (H) are 0.556 (T), 0.557 
(S), 0.563 (P) and 0.559 (H) respectively, (b) The slope (k) 
of line separating phase I and II versus the square root of the 
cross section area of nanodot with unit radius, (c) Phase dia- 
grams of cylindrical nanoparticles with different anisotropies. 
Solid squares represent cubic anisotropy (C) of different mag- 
nitude. Open circles with different colors represent uniaxial 
anisotropy (U) with Ka^ jD — 1 showing valid scaling behav- 
ior with r\ = 0.56. Solid triangles represent combination of 
both anisotropies (U-l-C) with Kc? jD = 1. (d) The slope k 
versus the strength of the uniaxial anisotropy. 

2(d)). Meanwhile, when both anisotropies are present, 
the slope is dominated by the uniaxial term. Hence an 
analysis of this slope can be used to determine the uniax- 
ial anisotropy experimentally, based on the information 
given in Fig. 2(d). 

Various lattice structures exists in nature. It is im- 
portant to know whether the scaling technique depends 
on lattice structure. We calculated the phase diagram 
for hep and face centered cubic (fee) lattices and their 
variance by rotating the lattice structure in the cylinder. 
The results remain invariant as long as all parameters 
{Ja? /D, Ka? /D and density of spins) are kept the same 
and X is not too small. The above results indicate that 
the scaling behavior is robust to details of lattice struc- 
ture, crystalline anisotropy, and geometric shape. 



B. Particles with Core Structure 

Interestingly, magnetic nanoparticles with core struc- 
ture exhibit complete similarity with respect to the lattice 
constant (See Fig. [3]). Similar effects have been reported 
by Landeros et. ali^ To analyze the effect of the core, 
we choose an ansatz {Sz = exp(— 2r^/?^)) introduced by 
Feldtkeller and Thomas^. We fit the results of Monte 
Carlo (MC) simulations with this ansatz and obtain ac- 
ceptable agreement (see the inset of Fig. 3). From di- 
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FIG. 3: (color online) Scaled phase diagram of a single- 
domain cylindrical magnetic nanoparticle taking the vortex 
core into consideration. The black hollow circles represent 
the phase diagram for cylindrical nanoparticle with core free 
model taken from Fig. 1. The scaling exponent rj = 0.5. The 
inset shows the fitting of the core function to the MC result 
for the case of J'/D = 100. 



mensional analysis, the core size 1//3 obeys a scaling law 
of the following form: 



1//3 = Lex<i>p{R/Lex,H/Lex,a/Lex), 



(8) 



where L^x is the magnetic exchange length, as before, 
and $/3 is a scaling function, independent of J, D. Nu- 
merically, we find that the scaling function $^ is approx- 
imately independent of all its arguments, giving roughly 
1//3 w O.QLex- We use this as an additional governing 
parameter in the numerical calculations. 

In the presence of the core, the critical height now 
satisfies a physical law of the form: 



H = giJ,D,R,a,l/f3). 



(9) 



From dimensional analysis, we find again only two inde- 
pendent governing parameters, define L^x, and write: 



H = L^x^g{R/Lex, a/ Lex, l/iPlex))- 



(10) 



Numerically, we find that $g approaches a constant as 
its second argument a/L^x becomes small. This is evi- 
denced by the collapse of the phase diagrams in Fig. 3 
with r] = 1/2. The collapse implies invariance under the 
transformation in Eq. O and thereby the complete simi- 
larity with respect to the lattice constant. However, this 
is consistent with the incomplete similarity with respect 
to an exhibited in the core-free approach, since we have 
an additional dimensional length 1//3 that plays the role 
of the lattice constant. Similar results obtain when we 
change the topology of the nanoparticle and introduce an 
inner radius. 

As to the phase diagram itself, the core stabilizes 
the vortex configuration significantly, pushing the phase 
boundary between FM and the vortex phase to smaller 



values of R and H by about 35%. 
affect Fig. 2(a)(c) as well. 



Similar effects would 



C. Cylindrical Nanorings 

Next we consider the effects of changes in topology on 
the phase diagram. More precisely, we investigate the 
phase diagram of hollow cylinders, i.e. nanoring struc- 
tures characterized by an inner radius Ri , an outer radius 
i?, and a height H. We find, as in the previous section, 
that the critical height exhibits complete similarity with 
respect to the lattice constant (77 w 1/2). This is a conse- 
quence of the additional length Ri that plays the role of 
the lattice constant in regulating the vortex core energy. 




FIG. 4: (color online) Phase diagrams of a cylindrical nanor- 
ing for two different x. There are two competing ferromag- 
netic phases at small {R, H) and a vortex phase at large 
[R, H). Because of the finite inner radius Ri, the onset of the 
phase transition line between the two ferromagnetic phases 
is shifted to finite values of R. Also, the vortex regime is 
more extended for larger Ri. The blue lines in the figure 
are guides to the eye, indicating that the triple points form 
approximately a straight line. 

Fig. [5] shows three-dimensional phase diagrams in the 
{Ri , i?, H) parameter manifold of the nanoring topology 
for two different values of the exchange couplings J'. 
Again, one observes two ferromagnetic regimes at small 
{R, H) values, competing with a vortex phase at larger 
{R,H). Moreover, one finds that for larger inner radii 
Ri the the vortex phase is more extended. This con- 
firms the idea that the ring structure stabilizes the vortex 
configuration. The reason for this is that the core area, 
which typically pays a high energy penalty, is deliber- 
ately avoided in the ring structure. Another new feature 
of these phase diagrams is that the line separating the 
two ferromagnetic phases is not straight anymore. In- 
stead, it now starts at finite R = Ri, and its slope changes 
smoothly to 1.81 as the ratio between R and Ri becomes 
very large. This relationship can be observed clearly in 
Fig. which is derived from Eq.(ll,13) in Ref. [H. 
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Here we calculate the relationship between the critical 
height Hc{R,Ri) as a function of (i? — Ri), resulting in 
the "star" symbols in Fig. EJb) which align exactly with 
the line of our numerical calculation. Finally, the most 
surprising feature of the phase diagrams in Fig. [4] is that 
the triple points {Rt,Ht) for different Ri approximately 
form a straight line indicated by the two blue lines. This 
property is shown more clearly in Fig. EJc), i.e. cylin- 
der height at the triple point (i?t) versus Ri. It gives 
us a critical Ric beyond which there exists no in-plane 
ferromagnetic phase. 
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FIG. 5: (color online) (a) For cylindrical nanorings, the phase 
transition line Hc{R,Ri) separating the two ferromagnetic 
regimes is not straight, in contrast to the topologically con- 
nected objects discussed above, (b) Phase diagram for a cylin- 
drical nanoring. {Ri = G.Snm, J/D = 5000). The data "A" 
represent analytical phase transition lines calculated from (a), 
which are observed to coincide with the numerical results, (c) 
Height at the triple point {Ht) versus inner radius (Ri). The 
best fit for x=0.1 is Ht = 16.5 - 5(±0.03) x Ri, whereas the 
best fit for x=0.06 is Ht = 12.7 - 4.95(±0.09) x Ri. Hence, 
the two lines are approximately parallel, and can thus be col- 
lapsed via scaling with rj — 0.51. (d) The triple point radius 
(Rt) versus x. 

Another observation worth mentioning is that the in- 
tersect phase diagram in the Ri — plane of the cylin- 
drical nanoring does not coincide exactly with the phase 
diagram of the simply connected cylinder (Fig. 1). It 
is closer to the case when the core structure is consid- 
ered (Fig. [3]). This phenomenon happens for the scaling 
exponent as well, which will be discussed right below. 

One last feature to be discussed here is that the line 
connecting Ht (see Fig. Uand Fig. [DJc)) is parallel for 
different values of the exchange coupling J'. Comparing 
the two phase diagrams for different J', we anticipate 
that there exists scaling behavior here as well, as long 
as all three coordinates {Ri,R,H) are scaled. However, 
some difficulties arise since Ri should be different for dif- 
ferent J's, meaning that one would need to know the 
scaling exponent rj in advance. Luckily, we can estimate 
the value of r] from Fig. \5lc) as the two straight lines 
should scale if there is a scaling behavior. Thus we first 



attempt to scale these two lines and find that they fit 
best when rj ~ 0.51. Then we use this rj to scale Ri and 
attempt to see whether the scaling behavior holds. Fig. 
[5jd) shows the result. The scaling exponent is 77 = 0.515 
which is within 1% of the estimated value 0.51. It is much 
closer to 0.5 in the finite core case, implying complete 
self-similarity, since an additional length Ri is added and 
neglecting the core structure in the vortex state has lit- 
tle effect for ring structure. With these results, we can 
easily calculate the critical inner radius Ric — llnm for 
the parameters we choose, above which a flat nanoring 
is always in the vortex phase. This is quite small com- 
pared to typical nanorings fabricated experimentall}»2^, 
and suggests that nanorings are generically in the vor- 
tex phase, since they are typically flat with height small 
compared with the width. 



D. Elliptically Shaped Particles 

It has recently been observed that there exists a double 
vortex conflguration in elliptically shaped ferromagnetic 
particles i2L22ii2£i^ The full phase diagram for this case 
as a function of height, semi-major axis {Ra) and semi- 
minor axis (i?b), however, has not yet been calculated. 
One of the difficulties to determine this phase diagram, 
using the technique outlined above, lies in finding an ad- 
equate paramctrization of the double vortex state. The 
naive approximation of two single vortices is far from sat- 
isfying (see Fig. 6(c)). As we will see below, the energy 
of two single vortices with discontinuous magnetization 
along the minor axis is significantly higher than of a true 
double vortex with continuously varying magnetization 
(see Fig. 6(a)). Without an accurate paramctrization of 
the double vortex one could only rely on Monte Carlo 
or micromagnctic simulations which are extremely time 
consuming, and this would make it impossible to obtain 
a complete phase diagram. 

Here we propose a simple function to parametrize the 
double vortex. In our Monte Carlo simulations, we ob- 
serve that the shape of the double vortex (Fig. 6(a)) 
looks much like the equipotential lines of two electric 
point charges with opposite signs placed at the centers 
of the vortex cores (Fig. 6(b)). By symmetry these cores 
should lie on the major axis of the ellipse. Let the dis- 
tance from the core centers to the center of the ellipse be 
Rq. Then the vector field S'(r) is given by 



where. 

Ex = 
Ej, = 



S{x,y) 



X — Rn 



-Eyi + Exj 
El+El 



(11) 



X -\- Ro 



[(X - i?o)2 y2]3/2 [(2. + jr?^)2 + y2]3/2 

y y 

[{X - RoY + y2]3/2 [(2; + Jl^Y + y2]3/2 
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FIG. 6: (color online) Double vortex configuration for J/D = 
10(a; = 0.002), Ra/Rb = 2 (arrows represent the directions of 
magnetization), (a) Monte Carlo simulation result, ec? jD = 
21.12 (b) Our parametrization, ea? jD — 21.11, .7^ between b 
and a is 0.990 (c) naive parametrization (two single vortices), 
ea^ / D = 20.91, T between c and a is 0.974. e is the energy 
per spin and JT is the fidelity defined in the text. 



Interestingly, the optimal positions of the vortex cores 
yielding the lowest energy configurations do not coincide 
with the ellipse foci, but are located at non-trivial po- 
sitions on the major axis with constant k = Ro/Ra- 
depends almost exclusively the aspect ratio (Ra/Rb), and 
depends only very weakly on size. Within the range we 
examined [H < AOnm,Ra < 30nm). k decreases by only 
2% as the size is increased. For different aspect ratios 
we find K = 0.44 ± 0.1. These values coincide with re- 
cent experimental results2Si21i^. We choose Ra/Rb = 2 
as an example. In this case, k = 0.44. To quantify 
the quality of our parametrization of the double vor- 
tex, we look at the energy per spin (e) and the fidelity 
= N~^^-Si ■ S'i, i.e. defined as the average dot 
product of spins on each lattice point of two configura- 
tions S{r) and S'{f). The energy of our parametrization 
(Fig. 6(b)) is significantly closer to the energy obtained 
by Monte Carlo (Fig. 6(a)) and its fidelity is significantly 
closer to 1 than the two single vortex parametrization 
(Fig. 6(c)). This is important because the energies of 
the single vortex and the double vortex configurations 
are very close. If one uses the naive parametrization, the 



double vortex could never be the ground state. 

Using the parametrization of the double vortex in 
EgfTTl we now apply the scaling procedure to obtain the 
phase diagram for elliptically shaped particles (sec Fig. 
7). Since there is no good description for the core of the 
double vortex yet, a core- free system is assumed for sim- 
plicity. We estimate that the boundary will shift to lower 
values of Ra and H by about 35% when taking the core 
into consideration. 




x=0.008 
x=0.010 
x=0.020 



200 300 

R (nm) 



400 



500 



FIG. 7: (color online) Scaled phase diagram of an elliptically 
shaped magnetic nanoparticle (Ka'^/D — 1 and Ja^/D = 
5000) as a function of its semi-major axis (-Ra) and height (H) 
with an aspect ratio 2. The four competing phases are (I) out- 
of-plane ferromagnetism, (II) in-plane ferromagnetism, (III) 
single vortex state and (IV) double vortex state. The scaling 
exponent is 77 = 0.55. 

As expected, the double vortex state becomes sta- 
ble when both the semi-major axis and height of the 
nanoparticle are increased. In the vicinity of the phase 
boundary between the single vortex and the double vor- 
tex states, the energies for the two configurations are very 
close, and hence there could be a large metastable region 
close this phase boundary where both states could exist 
in nature. This is likely the reason why both these config- 
urations have been observed in experiments on the same 
particle^^. Regarding the scaling exponent, 77 = 0.55 is 
again observed in this core-free consideration, implying 
incomplete self-similarity. 

Here we have only focused on the double vortex state. 
When the system size and the aspect ratio are sufficiently 
large, it is possible that multivortex states emerge. Be- 
sides such complex single domain structures, cross-tie do- 
main walls^ could exist in these structures as well. It 
would be highly interesting to know under which condi- 
tion these configurations could be stabilized. 



IV. CONCLUSIONS 

In conclusion, we have extended and analyzed the hy- 
pothesis of physical similarity put forward in Rell6l. Re- 
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gardless of shape, anisotropy, or crystal structure, wc find 
numerical evidence for incomplete similarity (77 = 0.55) 
with respect to the lattice constant a when a "core-free" 
model is assumed. Introducing additional small length 
scales, such as core-size or an inner radius, restores com- 
plete similarity (7; ~ 0.5), since the new small length 
regulates the vortex core. 

A three-dimensional phase diagram for the cylindrical 
ring structure was obtained and a linear relationship be- 
tween the height (-ffj) at the triple point and the inner 
radius (Ri) was found, which offers a straightforward way 
to calculate the critical inner radius above which there ex- 
ists no in-plane ferromagnetic phase. A new parametriza- 
tion for double vortex configurations was proposed. This 



configuration was found to be the ground state when both 
the radius and height of the cUiptically shaped magnetic 
particle are large. Finally, a new phase diagram for ellip- 
tical nanoparticles including a double vortex phase was 
determined. 
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